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Within a Monte Carlo simulation we investigate the statistical properties of an electron flow 
injected with a Poissonian distribution and transmitted under ballistic regime in the presence of 
long-range Coulomb interaction. Electrons are shown to exhibit a motional squeezing which tends to 
space them more regularly rather than strictly at random, and evidence a sub-Poissonian statistics 
with a substantially reduced Fano factor F n <C 1. The temporal (anti)correlation among carriers is 
demonstrated to be a collective effect which persists over the transit of several successive electrons 
and results in a considerable (more than one order of magnitude) shot noise suppression. 



The study of electronic noise in low-dimensional sys- 
tems, and mesoscopic conductors in particular, Jaas been 
receiving increasing attention in recent years.tFLl Indeed, 
by complementing the analysis of current- voltage charac- 
teristics, the knowledge of statistical properties provides 
valuable information on the microscopic mechanisms ul- 
timately responsible for charge transport. Uncorrelated 
carriers exhibiting Poissonian statistics are known to be 
characterized by a shot noise power given by the Schot- 
tky formula Si = 2ql (q being the unit charge and / the 
average current). To a general extent the ideal 2ql value 
can be reduced by imparting an anticorrelation to the 
electron stream, thus making the electron statistics be- 
come sub-Poissonian. This can be achieved by exploiting 
several mechanisms. Pauli correlation under degenerate 
conditions has been investigated in recent experiments!!] u 
and theories (see detailed bibliographyO). A 1/2 noise- 
suppression factor |has been found in a symmetric double- 
barrier junction.ldi3 For a diffusi-ve conductor a 1/3 sup- 
pression has been demonstrated!] The statistics of charge 
transmitted through _a degenerate conductor has been 
found to be binomial,!! while that through a double bar- 
rier junction., is in between a Gaussian and a Poisson 
distribution.^ 

Another possibility to incorporate correlations among 
electrons is by means of long-range Coulomb repulsion. 
The desired effect can be achieved by virtue of a poten- 
tial barrier near an injecting contact, which fluctuates 
synchronously with electron passages through it. In this 
way, an incoming Poissonian flow is converted into an 
outcoming sub-Poissonian flow. This effect is similar to 
that leading to shot-noise suppressiop-in vacuum diodes 
and discussed in literature long agojld although, in the 
case of a device, the potential barrier should not neces- 
sarily be provided by space charge. In a semiconductor 
structure it might be, e.g., a constriction, as in the ex=. 
periment on a quantum point-contact by Reznikov et alu 
We notice that although the importance of Coulomb cor- 



relations on noise suppression in mesoscopic structures 
has been.^pphasized,Ei}E3 only a few calculations are 
available, E£lEil and the statistics of charge transmitted un- 
der Coulomb correlations has never been studied to the 
best of our knowledge. 

It is the aim of the present paper to address the prob- 
lem of Coulomb correlation by investigating the carrier 
statistics under nondegenerate ballistic transport regime 
through a Monte Carlo simulation. We have found the in- 
teresting result that just after the potential barrier, self- 
consistently generated by the electron flow, the statis- 
tical variance of electron-number counts decreases with 
the distance from the barrier. Thus, correlated electrons 
exhibit a motional electron-number squeezing, character- 
ized by a reduced Fano factor F re < 1. 

For our analysis we consider the following simple 
model: a lightly doped active region of a semiconductor 
sample sandwiched between two heavily doped contacts 
(of the same semiconductor) subject to an applied volt- 
age bias U. The contacts are considered to be ohmic 
(the voltage drop inside them is negligible) and they re- 
main always at thermal equilibrium. Thus, electrons are 
injected from the two opposite contacts according to a 
Maxwell-Boltzmann distribution at the lattice tempera- 
ture Ti,. To exclude correlations due to Fermi statistics, 
the electron gas is assumed to be nondegenerate. Elec- 
trons are injected with a Poissonian statistics, i.e., the 
time between two consecutive electron injections is gen- 
erated with a probability per unit time P(t) — Toe -1 " ', 
where Tq = ^n c vthS is the injection-rate density, n c the 
carrier conce ntration at th e contact, S the cross-sectional 
area, i>th = y // 2ksTi J /Trm the thermal velocity, and m the 
electron effective mass. For the simulation we used the 
following set of parameters: Tl=300 K, to = 0.25 too 
(too being the free electron mass), relative dielectric con- 
stant e=11.7, active region length L = 200 nm, and 
n c = 4 x 10 17 cm~ 3 , much higher than the sample dop- 
ing (here taken of 10 11 cm~ 3 , but the same results are 



obtained up to 10 15 cm -3 ). Under these conditions the 
electron dynamics is determined by the ratio A = L/Lu c , 
where Lr> c is the Debye screening length associated with 
the electron injection density. The above set of parame- 
ters corresponds to A = 30.9. 

We calculate the fluctuating current I(t) and its low- 
frequency fluctuation spectrum Si by an ensemble Monte 
Carlo simulator self-consistently coupled with a Poisson 
solver (PS).EJ To illustrate the importance of Coulomb 
correlation we provide the results for two different simu- 
lation schemes. The first one involves a dynamic PS: the 
fluctuations of the potential follow self-consistently the 
fluctuations of the space charge by virtue of the Pois- 
son equation. The second scheme makes use of a static 
PS, which calculates the stationary potential profile only, 
so that carriers move in the frozen non-fluctuating elec- 
tric field. Both schemes give exactly the same average 
current and steady-state spatial distributions of all first- 
order quantities, but they provide different noise levels 
and statistical distributions of transmitted carriers. 

Fig. |] illustrates the potential minimum which, act- 
ing as a barrier for the electrons moving between the 
contacts, is shifted by changing the bias U. The cor- 
responding noise-suppression factor 7 = Si/2ql is re- 
duced considerably when the dynamic PS is used, while 
for the static PS 7 w 1 for all the biases (see Inset of 
Fig. ^) . The noise suppression becomes particularly pro- 
nounced at qU « 40fcsTi, when the potential barrier al- 
most vanishes completely and its fluctuations affect the 
transmission of the most populated (low-velocity) states 
of injected carriers. At the highest voltages the bar- 
rier disappears and no shot-noise suppression is observed. 
Note that the shot-noise suppression factor by Coulomb 
interaction can achieve as low value as desired in coru 
trast to other mechanisms (1/3 £pr diffusive conductorsjD 
1/2 for double-barrier junctionsu) by increasing the ratio 
A = L/Ld c , i.e., by increasing the sample length and/or 
the carrier concentration at the contact, provided the 
transport remains ballistic. In the following we fix the 
bias at U = 40k B T L /q, for which 7 w 0.045. 

Having found the level of noise below the ideal Pois- 
sonian value, it is natural to ask which are the statisti- 
cal properties of the transmitted carriers once they are 
temporally correlated. To address this question we reg- 
ister the times of passage of electrons through different 
sections of the sample at a distance x from the injecting 
contact. We register only those going from the cathode to 
the anode, since these are the only electrons contributing 
to the low-frequency current noise. With this procedure 
we are able to study how the carrier statistics imposed 
at x=0 (in our case Poissonian) modifies at different dis- 
tances from the cathode. For each section x the times of 
passage of carriers {t±, £2, • • • , U, ■ ■ •} constitute a station- 
ary stochastic point process, which can be described by 
specifying the random set {tj}, where r< = —t% is the 
time interval between events. By knowing the set {t;}, 
we calculate the distribution function P n (T), which is the 
probability of detecting n electrons during the observa- 
tion time interval T . For a Poissonian process all time 
occurrences are statistically independent, which leads to 



the simple formula 

Pn(T) = 5^e- rT , n = 0,1,2,... (1) 

where V is the rate density of events. The distribution (|l|) 
is characteristic of uncorrelated transport and is tested 
to correctly describe the carrier statistics at the injecting 
contact at x=0. The distribution of carriers at the re- 
ceiving contact depends crucially on the PS scheme (see 
Fig. |). For the static PS P n {T) is perfectly fitted by the 
same Poissonian formula (|l|) as for the injected carriers, 
but with a smaller value of the rate density Y = kYq, 
with k ~ 0.65. Its reduction is caused by the part of 
injected carriers which is-^eflected by the potential bar- 
rier back to the contact.Ej In contrast, for the dynamic 
PS the distribution function no longer obeys the formula 
(|l|). Fig. [2] shows the following differences between the 
dynamic and static PS: (i) for each n > 1 the maximum 
of the distribution is shifted to longer T; (ii) the proba- 
bility distribution is narrowed; (iii) the higher the index 
n, the more the dynamic case deviates from the static 
one, the distribution profiles becoming more symmetri- 
cal and closer to a Gaussian shape. We remark that the 
difference at point (ii) can be interpreted as a motional 
squeezing of electron number and corresponds to a higher 
regularization of the carrier passage due to correlation 
among them. 

The temporal correlation among carriers can be char- 
acterized by the autocorrelation function of the inter- 
event times as a function of the number of successive 
events defined as 

C T (k) = ((r j -r)(r j+k -T)), fc = 0,l,... (2) 

Here r is the mean of the random variable Tj, and an- 
gular brackets denote average over the reference index j. 
Fig. H shows the function C T {k) calculated for the static 
and dynamic PS using the random sets {r^} registered at 
different positions x along the sample. Firstly we notice 
that for the static PS the correlation function is zero for 
all k 7^ [Fig. ^(b)]. This means that the absence of 
correlation among carriers imposed by the Poissonian in- 
jection is preserved at any section of the sample. Indeed, 
being the electric potential frozen no additional correla- 
tion is introduced inside the system. The only nonzero 
value is the variance C r (0) = Var(r) = r 2 — r 2 , which 
increases with x until x coincides with the minimum of 
the potential profile (here the carrier flow is decreased 
due to the reflection), and then saturates for larger val- 
ues of x (see Fig. [|). In contrast, for the dynamic PS 
C T (k) is in general nonzero for k ^ and, moreover, it 
is negative by decaying to zero at increasing values of 
the index k [Fig. ||(a)]. Negative values of C T (k) start 
to appear just after the potential minimum, their magni- 
tude tending to increase for longer distances. The nega- 
tive sign means that once a time interval r m longer than 
the average value r is registered, to compensate such a 
deviation the successive r m+1 is expected to be shorter 
than t. The mechanism responsible for this compensat- 
ing behavior is the self-consistent redistribution of the 
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potential which, by varying the barrier height introduces 
a correlation into the carrier flow. From Fig. |^(a) we 
estimate that the number of consecutive carrier passages 
over which the correlation persists is ~ 6 [i.e., C T (k) p» 0, 
for k > 6]. It should be stressed that the mean time 
between consecutive electrons is not influenced by the 
long-range Coulomb correlation, since r is related to the 
average current. Fig. || shows that for both the static 
and dynamic PS r increases in the region before the bar- 
rier (x < 0.06L), due to the loss of the part of injected 
electrons which are reflected by the potential barrier, and 
then saturates. The Coulomb correlation is mirrored by 
the second moment Var(r) which, for the dynamic PS, 
is found to decrease with distance just after the barrier 
(see Fig. f|). 

The autocorrelation function C T (k) refers to the statis- 
tical distribution of inter-event times {r^} which, in turn, 
is related to the operation mode of registering the carrier 
arrival times. Another operation mode is the number- 
counting, for which the number of carriers nf in a given 
time interval T is registered. For this mode the statisti- 
cal information is carried by the random variable {nf} 
which depends on the observation time T as a parame- 
ter. The first-order statistics of the random set {nf} is 
described by the distribution function P n {T) shown, for 
example, in Fig. || for the carriers at the receiving con- 
tact. A better insight into the carrier correlation can be 
obtained by calculating the correlation function 



C n (k,T) 
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n T )}, k = 0,l,..., (3) 



representing the second-order statistics. Here, n T is the 
mean number of counts during the observation time inter- 
val T. For fixed T ~ r, C n (k, T) behaves like C T (k), i.e., 
it is negative for k > 1 with its lowest value at fe=l [like 
in Fig. 0(a)]. Therefore, C„(l, T) can be a good measure 
of the correlation between the number of counts. This 
quantity, normalized to the mean number of counts n T , 
is shown in Fig. || for the statistics registered at the re- 
ceiving contact under dynamic regime. C„(1,T) is seen 
to display a minimum as a function of T (indicating the 
maximum of correlation) for T « 5/r. This means that 
the maximum of the correlation occurs when the obser- 
vation time interval T is chosen such that the number of 
counts is n T w 5. In fact, this is another estimation of 
the number of consecutive carriers over which the cor- 
relation persists, in addition to that obtained from the 
approximate decaying range ofCV(fc) ofFig.|(a) (A; « 6). 

The correlation function for k=0 gives the variance of 
the number of counts C„(0,T) = Var(n T ) = (n T ) 2 — 
(n T ) 2 . In Fig. || we present also the ratio of the vari- 
ance to the mean F n (T) — V&r(n T )/n T kaown as the 
Fano factor in photon counting statistics.ll3 It charac- 
terizes the degree of deviation of the counting statistics 
from the Poissonian one, whenever F n (T) ^ 1. We have 
found that for the static PS F n (T) = 1 always, while for 
the dynamic PS it decreases with T (see Fig. ||). Notice 
that if the observation time T < 1/T, F n (T) w 1, i.e. 
Coulomb correlation among carriers cannot be detected 
by any means, and the observed statistics of carriers is 



always Poissonian. For T ^> 1 /V the correlation between 
the counting intervals is lost [C„(1,T) — » 0], while the 
Fano factor approaches its minimum value at the shot- 
noise suppression factor F n (oo) = 7 w 0.045 (Fig. |). 

In conclusion, within a Monte Carlo scheme we have 
investigated the statistics of carriers which are Poissonian 
injected and transmitted through a nondegenerate ballis- 
tic structure in the presence of long-range Coulomb inter- 
action. The sub-Poissonian statistics exhibited by carrier 
number fluctuations registered at the receiving contact 
provides an interesting example of motional squeezing in 
analogy with photon satististics. The negative value of 
the correlation function of both: time intervals between 
consecutive electrons passage, and number of counts in 
successive observation time intervals, mirrors the tempo- 
ral correlation of transmitted carriers introduced by the 
presence of space charge. The small value of the Fano 
factor (F n <C 1) here found, besides being of fundamen- 
tal interest, offers interesting perspectives in communica- 
tion systems. Indeed, a sub-Poissonian electron flux can 
generate a photon-number-squeezed light which, in turn, 
can be used to enhance a channel capacity to-transmit 
information in optical communication systems.113 
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FIG. 1. Spatial profiles of the normalized potential (f)(x) 
for different applied voltages U (in units of fcsTi/g). Inset: 
shot-noise suppression factor vs voltage U for the dynamic PS 
(7dyn) and static PS (7st), respectively. 
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FIG. 2. Distribution function P n (T) of electron counting 
statistics at the receiving contact for the static PS (dotted 
line, Poissonian statistics) and for the dynamic PS (solid line, 
sub-Poissonian statistics) . 




FIG. 3. Correlation function of inter-event times as a func- 
tion of the number of successive events at different distances 
x from the injecting contact for the dynamic (a) and static 
(b) PS, respectively. For a better view discrete points are 
connected by lines, x is in units of the sample length L. 
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FIG. 4. Mean and variance of inter-event times as a func- 
tion of distance x from the injecting contact, r is multiplied 
by the injection rate To. 
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FIG. 5. Fano factor F n (T) (left axis, squares), and correla- 
tion function of the number of counts between two consecutive 
observation intervals, C„(l,T), normalized by n T (right axis, 
stars), calculated at the output contact within the dynamic 
scheme. 
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